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Context. Massive, eclipsing, double-lined, spectroscopic binaries are not common but are necessary to understand the evolution of massive 
stars as they are the only direct way to determine stellar masses. They are also the progenitors of energetic phenomena such as X-ray binaries 
and y-ray bursts. 

Aims. We present a photometric and spectroscopic analysis of the candidate binary system Cyg OB2-B17 to show that it is indeed a massive 
evolved binary. 

Methods. We utilise V band and white-light photometry to obtain a light curve and period of the system, and spectra at different resolutions to 
calculate preliminary orbital parameters and spectral classes for the components. 

Results. Our results suggest that B17 is an eclipsing, double-lined, spectroscopic binary with a period of 4.0217 + 0.0004 days, with two 
massive evolved components with preliminary classifications of 07 and 09 supergiants. The radial velocity and light curves are consistent with 
a massive binary containing components with similar luminosities, and in turn with the preliminary spectral types and age of the association. 



Key words, binaries: eclipsing - 
Cyg OB2B17 



1. Introduction 



binaries: spectroscopic - stars: early-type - stars: fundamental parameters - stars:evolution - stars: individual: 



O stars are amongst the most massive and intrinsically lumi- 
nous stellar objects found in galaxies. Since they are the only 
direct way to measure the masses and radii of stars, binaries 
are the perfect testbeds for studying the physical prop erties and 
evolut ion of such stars. Unfortunately, as reported bv lBonanos 
d2009h . less than 20 O stars have accurate (<10%) dynami- 
cal mass estimates. Because of the scarcity of known massive, 
eclipsing, double-lined, spectroscopic binaries - and hence dy- 
nami cal mass e stimates for stars at different evolutionary states 
(e.g., Gies 2002) - the mass luminosity relation and theoretical 
evolutionary tracks of massive stars (M > 2OM ) are currently 
poorly constrained by observations. 



In order to address this shortfall, much effort has been ex- 
pended to identify further examples, utilising photometric and 
specnosco pic observations of y oung massive clusters such as 
the Arches dMartins et alJ2008h. Quintuplet dFiger et al.l 1999b . 
and Westerlund 1 ddark et al.l2005HRitchie et al.l2009h . These 
observations indicate that the binary fraction is potentially very 



higrQ; iKobulnickv & Frvetl d2007t) i nferred it to be > 80% 

~ d2008h esti- 



in Cygnus OB2 (Cyg OB2), while Clar k et al. 



mate the binary fraction of WR stars in Westerlund 1 to be 
>70%. Currently the stars with the highest dynamical mass 
estimates are the twin WN6 ha components of WR20a within 
Weste rlund 2 (83M m +82M m : iRauw et alJl2004 iBonanos et al ' 



|2004 and the newly disc overed WN6ha bina ry NGC3603-A1 
(116+31M Q + 89+16MM: ISchnurr et alj2008h . This in turn has 
implications for the determination of the Initial Mass Function 
(IMF) for the clusters in question and, by extension, the e mpir- 
ically determined maximum mass possible for a star (cf. Figer 



2005). Moreover, massive close binaries are the progenitors of 



such diverse energe tic phenome na as supernovae, y-ray bursts 
and X-ray binaries (Ribas 2006). Clearly the properties of the 
progenitor binary population must be known to constrain their 
formation channels. 

Cyg OB2 is one of the most massive and richest asso- 
ciations in the Gala xy. It is ~2 Myr old and 1.8 kpc away 
(Ki minki et al] 120071) . Containing at least 60-70 O-type stars 



1 A large amount of observations and patience are necessary to de- 
termine the true binarity of a sample of stars in an open cluster as 
observed in lSana et all b008h 
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(Neg ueruela et al.ll2008l) . its proximity and accessibility to op- 
tical studies have made it the focus of numerous observational 
campaigns to d etermine the properties of its massive stellar 
population (e.g. Massey & Thomson 1991; Knodlseder 2000; 
Hansod2003HComer6n et al.l2002tlKiminki et al.l200m2008). 



Cyg OB2 B17 dComeron et al.1 120021: henceforth B 17 and 
also known as V1827Cyg, 2MASS J20302730+41 13253, 
NSVS 5738756; a 2 ooo=20 /, 30'"27.3 J , 5 2 ooo=+41"13'25"; 
V=12.6) is a luminous, variable member of the Cyg OB2 

(2002) observed the system as 



Comeron et al 



association. 

part of a near-infrared spectroscopic survey and found the 
spectrum presented Bry emission, confirming it to be an 
evolved massive star. F ollow up observations were made by 
Negueruela et al.l (2008), who found it to stand out from the 
rest of the members due to its variability and strong emission 
lines. They classified it as an Ofpe star and suggested it was a 
strong binary candidate. 

This paper reports the first results of an extensive multi- 
epoch photometric and spectroscopic observational cam- 
paign on the binary candidate B17. Section 2 gives a de- 
scription of the photometric and spectroscopic observations. 
Photometrically, the system was found to be variable and we 
report the analysis of the light curve in Section 3. More spectro- 
scopic data were obtained permitting preliminary spectral and 
luminosity classification of the system; this analysis is shown 
in Section 4, along with descriptions of the long and short 
timescale variations of the spectra. The light and radial velocity 
curve modelling is presented in Section 5. A discussion of the 
system, including its evolutionary status, is found in Section 6 
and a summary is presented in Section 7. 

Note that the central goals of this manuscript are to verify 
the binary hypothesis and present a preliminary spectral clas- 
sification. A full analysis of the system, consisting of the de- 
convolution of an expanded spectral data set and subsequent 
model atmosphere analysis to determine the fundamental stel- 
lar parameters of the system will be presented in a future paper 
(Stroud et al. in prep). 



2. Observations and Data Reduction 

2.1. Photometry 



The North Sky Variability Survey (NSVS; IWozniak et al 



2004) is a record of the sky at declinations higher than 5 
-38° over the optical magnitude range 8 to 15.5. It contains 
light curves of over 14 million objects. The data were taken 
between 1999 April and 2000 March by the first-generation 
Robotic Optical Transient Search Experiment (ROTSE-I) at 
Los Alamos National Observatory, New Mexico. The telescope 
consisted of four unfiltered Canon 200mm telephoto lenses 
with f/1.8 focal ratio, each covering 8^2 x 8*?2. These were 
equipped with APlO-cameras and Thomson TH7899M CCDs. 
The lenses had a typical point-spread function with a full width 
half maximum of ~20". In a median field, the bright unsat- 
urated stars had a point-to-point photometric scatter of -0.02 
mag and position errors within 2". T he calibrated images wer e 
passed through SExtractor software dBertin & Arnouts 19961) . 
reducing them to object lists. The data were accessed through 



the Sky Database for Objects in Time-Domain (SkyDOT) at 
Los Alamos National Laboratory. A total of 186 observations 
were obtained for B 17. 

Additional V band photometry was obtained by amateur as- 
tronomers Pedro Pastor Seva (observer 1) and Manuel Mendez 
Marmolejo (observer 2) between 2007 April 03 and 2007 June 
11. Observatory 1 is located in Muchamiel (Alicante, Spain). 
The telescope used was an 8 inch Vixen VISAC Schmidt- 
Cassegrain telescope. It has a field of view of 24'x20' and 
a focal ratio of f/8.The telescope was equipped with a SBIG 
ST10-XME CCD chip. Observatory 2 is located in Rota (Cadiz, 
Spain). The telescope used was an 8 inch Meade LX200 
Schmidt-Cassegrain telescope. It has a field of view of 16'xl2' 
(as a smaller SBIG ST7-XME chip was used) and a focal ra- 
tio of f/6.3. Exposures varied between 3 and 4 minutes giv- 
ing a SNR ~200. The data were calibrated and analysed using 
Mira Prc@ and AIP4WINv^| packages. The apparent magni- 
tudes were obtained using differential photometry with respect 
to a set of reference stars in the image with known magnitudes. 
In both cases the precission for the individual meassures is al- 
ways better than 0.01 mag. 



2.2. Spectroscopy 

The spectra were obtained from several telescopes during the 
course of ~three years. Table Q] lists the full set of observa- 
tions. The first set was observed with the 1.52-m G. D. Cassini 
telescope at the Loiano Observatory (Italy) during the night of 
2004 July 18. The telescope was equipped with the Bologna 
Faint Object Spectrograph and Camera (BFOSC) and an EEV 
camera. Grism 3 was used, which covers 3300-5800A with 
a resolution of ~6A. More spectra of the system were ob- 
tained with the 4.2-m William Herschel Telescope (WHT), in 
La Palma (Spain) equipped with the ISIS double-beam spectro- 
graph, during service runs on 2006 June 1 1 and August 18. The 
instrument was fitted with the R300R grating and MARCONI2 
CCD in the red arm and the R300B grating and EEV 12 CCD 
in the blue arm. Both configurations result in a nominal disper- 
sion of 0.85 A/pixel (the resolution element is approximately 3 
pixels in the blue and 2 pixels in the red). 

The system was also observed with the 2.5m Isaac Newton 
Telescope (INT) in La Palma on 2006 September 10-11. 
The telescope was equipped with the Intermediate Dispersion 
Spectrograph (IDS), fitted with a R632V grating and EEV10 
CCD. Another 12 spectra were obtained during a dedicated 
run on 2007 August 21-22 at the WHT. The system was ob- 
served in the blue arm with grating R1200B (nominal dis- 
persion of ~0.23A/pixel). All the spectra were re duced with 
the Starlink packages c cdpack dDraper et al.l 120001) and figaro 
( Shortridee et al.l 19 97) and analysed using figaro and dipso 
(Howarfh et al.lll998l) . 



2 Mira Pro software is published by Mirametrics Inc., which has no 
connection with the Monterey Institute for Research Astronomy 

3 AIP4WINv2 is published by Willmann-Bell, Inc. 
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Table 1. Log of Spectroscopic Observations. 



Date 


Telescope 


Nominal 


Wavelength 


Phase 




/Instrument 


Dispersion 
(A/pix) 


Range (A) 




04/07/18 


Cassini/BFOSC 


6 


3800-6400 


0.437 


04/07/18 


Cassini/BFOSC 


6 


6100-8200 


0.450 


06/06/11 


WHT/ISISB 


0.86 


3200-5200 


0.036 


06/06/11 


WHT/ISISR 


0.93 


5400-8000 


0.036 


06/08/18 


WHT/ISISB 


0.86 


3200-5200 


0.929 


06/08/18 


WHT/ISISR 


0.93 


5400-8100 


0.931 


06/09/10 


INT/IDS 


0.9 


3900-5000 


0.600 


06/09/11 


INT/IDS 


0.9 


3900-5000 


0.840 


06/09/11 


INT/IDS 


0.9 


3900-5000 


0.891 


07/07/21 


WHT/ISISB 


0.86 


3600-5300 


0.738 


07/07/21 


WHT/ISISR 


0.93 


5400-8300 


0.738 


07/07/21 


WHT/ISISR 


0.93 


5400-8300 


0.740 


07/07/21 


WHT/ISISR 


0.93 


5400-8300 


0.741 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.380 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.385 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.397 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.416 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.436 


07/08/21 


WHT/ISISB 


0.23 


3900-4750 


0.453 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.629 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.634 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.648 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.664 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.683 


07/08/22 


WHT/ISISB 


0.23 


3900-4750 


0.701 



3. Photometry 

The Starlink software period ( Dhillon et al.ll200 1 ) was used on 
both the NSVS and amateur photometry to search for a mod- 
ulation period in the photometric data, using phase dispersion 
minimisation, x 2 of sine fit vs frequency and string-length vs 
frequency methods. The resultant periods were consistent and 
favoured an orbital period of 4.0217 + 0.0004 days (the error 
being estimated from the spread in the periods derived from the 
different methods employed). When both data sets were folded 
together they were found to be in phase; therefore, given that 
there appears to be no change of shape or shift in period in the 
7 years between the NSVS and amateur observations, we are 
confident that the period determined is accurate to within the 
errors quoted. 

Fig. Q] shows the periodogram obtained using Period with 
the reduced-^ 2 technique (top) and the V band light curve 
folded on a 4.0217 day period, along with the phases of the 
spectra obtained (bottom). 

The shape of the light curve suggests that the system is a 
semi-detached binary; Both minima are narrow and demon- 
strate different eclipse depths, although the 0.5 - 0.6 mag range 
indicates the almost complete eclipse of each star and hence 
that the two stars are of similar size but with different luminosi- 
ties (in a contact system, the temperature of both stars should 
be the same) There appears to be an asymmetry on the light 
curve betwe en phases 4> = 0.6 - 0.75 which is observed in both 



o 



0.2 0.3 

Frequency (1/days) 




Fig. 1. Upper panel: Periodogram for B17 data using the 
reduced-^ 2 technique. Lower panel: Light curve folded on a 
4.0217 day period. The red triangles show the phases for which 
spectra were obtained. 



(2009) round a similar depression in the light curve for the bi- 
nary LMC-SC1-105. They attributed t hese asymmetries to the 
presence of a mass-transfer stream. Linde r et al. I d2009h also 
observed this asymmetry in the overcontact binary Cyg OB2 
#5 and attributed this to the secondary being brighter on its 
leading side due to the colliding wind region; first discussed by 
Rauwet al.l dl999t) . 

The quadratures display an O'Connell effect - where the 
maxima are of different brightnesses - of ~0.02 mag which is 
likely due to variations in the brightness of the stellar surface(s) 
due to mass transfer between the components. 

The principal minimum was estimated using the amateur 
data to be at JD = 2,454,272.527 + 0.005. From the light 
curve, the ephemeris for the primary eclipse was found to be: 



mini = 2454272.527 + 4.02 17 E (JD) 



(1) 



sets of data. Hildit ch et al.l ([2005) found similar depr essions in 
eclipsing binaries in the Small Magellanic Cloud and lBonanos 



where E is the number of orbital cycles after the given epoch. 

We note that our period dete rmination is in agreement with 
a value recently established by Oterol ( 20081) using the same 
data from NSVS as used in this paper. 
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4. Spectroscopy 

The spectrum of B17 is found to be highly variable on both 
short (<day) and long (~year) timescales, with the most promi- 
nent features in the spectra being H, He i, He n and N m lines 
in absorption and He n and N m lines in emission. Wavelength 
shifts are observed for different lines, some of which become 
double in some spectra suggesting that it is a binary system 
with at least one hot component, as revealed by the presence of 
He ii and N m. 



4.1. Long Term Variability 

Fig. 2 shows blue-violet spectra (4000-4900A) obtained at ran- 
dom phases between 2004 July and 2007 August. The most 
prominent elements observed are H, He i, He n and N m in ab- 
sorption and He n ,14686 and N m ,1,14634-40-42 in emission. 
The spectra show strong, highly variable N in and He n emis- 
sion lines; in 2004 July 18 these are of a similar strength but by 
2006 June 1 1 the He n line is almost three times the intensity 
of the Nm lines. The Nm profile is also variable; it appears as 
two separate peaks in some of the spectra (2006 June 11, 2006 
September 10-11, 2007 July 21) and it is almost completely 
blended in the spectra taken on 2004 July 1 8 and 2006 August 
18, although this in part could be due to the low resolution of 
the spectrum from 2004 July 18. 

The hydrogen lines are also highly variable; this is partic- 
ularly evident for Hy when compared to the diffuse interstellar 
band (DIB) at ,14428. H/3 shows a P-Cygni profile varying in 
both width and strength suggesting the stellar wind is highly 
variable. Finally, the Hei ,14471/ Hen ,14541 ratio, which is 
empl oyed as a standard diagnostic f or the spectral type of early 
stars (IWalborn & Fitzpatricklll990h . is also observed to vary 
over the course of the observations. 

The most prominent line in the yellow spectra (5600- 
6000A) is the varying P-Cygni profile in Hei /15875, which 
blends with the Nai interstellar lines (Fig. 3). When compared 
to the nearby DIBs at 5780-5800A, the absorption trough in 
He i line was observed to vary by a factor of 4 in intensity rel- 
ative to the nearby DIB features, with the emission component 
also varying in strength by a factor of 2. Other lines in the yel- 
low spectra include C m ,15696 and C iv around ,158 10. 

The most prominent feature in the red spectra (6350- 
6650A) is the Ha line which also demonstrates a highly vari- 
able P-Cygni profile which is blended with He n absorption fea- 
tures in all the red spectra, complicating analysis of the profile. 
It has a double peak profile which varies in intensities with 
the redder peak being between 2 and 4 times stronger than 
the bluer peak. The depth of the He n line, which presumably 
contributes to the double peaked morphology is also variable 
(Fig. 4). Preliminary calculations for the velocities of the blue 
edge of the Ha profile (v e dge) have a range of 2260 - 2625 
kms~' with an average 2440+50 kms -1 . The v e d ge is related 
to the terminal velocity of the wind; in the extreme ultravio- 
let, the terminal velocities for OB star s are 15%-20% smaller 
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than the edge velocities (Pri nia et al 



wind profiles of OB stars are thought to be a consequence of 
highly structured winds on multiple scales, and while it appears 



Fig. 2. Blue-violet spectra of B17 obtained at different phases 
over a three year period, with the most prominent lines labeled. 



likely that the variations observed for B 17 arise due to a highly 
anisotropic circumstellar envelope observed at differing lines of 
sight throughout the orbital period, the current limited data set 
offers little prospect of a more explicit physical interpretation. 

Unfortunately, no spectra with the same orbital phases have 
been obtained at different epochs and so it has not been possible 
to search for unambiguous long term secular variability in the 
current data set. 

In Fig. 5 we plot the EW of the He n ,14686 line against 
orbital phase, finding the EW to be broadly anticorrelated with 
the photometric light curve, increasing by a factor of two dur- 
ing the entry to eclipse. This result suggests that the strength 
of the line remained -constant over the segments of the orbital 
period sampled by the observations, implying that the change 
in EW is primarily due to dilution by the variable continuum. 

4.2. Short Term Variability 

The set of 12 high resolution spectra obtained during the nights 
of the 2007 August 21-22 cover the phases 0=0.38-0.453 and 
0=0.629-0.701 of one period. The spectra obtained during the 
first night (Fig. 6 upper panel) have narrower H and He n ab- 
sorption lines lines compared to the spectra from the second 
night. During the first night, the Hei ,14471 /Hen ,14541 ratio 
- used for classifying O stars - changes monotonically in the 
sense of increasing He i ,1447 1 strength, reaching unity in the 
final spectrum. 
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Fig. 3. B17 yellow spectra showing variability at different 
phases taken over a period of three years. Prominent lines are 
labeled. 

A second, blueshifted absorption component appears on the 
shoulder of the He i A4471 line, becoming stronger with time; 
similar evolution is also observed in the US profile. Initially, 
the N m /14634-40-42 emission is highly blended, with the red- 
shifted peak being the strongest. These peaks start to separate 
on the last spectrum of that night (0=0.453). The Hen /14686 
emission feature appears to be non-symmetric and strengthens 
throughout the night. The He i i4026, He n /I4200 and Hy lines 
do not vary significantly throughout the first night. 

During the 17hrs between the last observation of night one 
(0=0.453) and the first of night two (0=0.629) the spectrum has 
clearly evolved. The He i i4026, He n /14200 and Hy lines are 
broader and non symmetrical. On the second night (Fig. 6 lower 
panel), the He i A4471 / He n A4541 ratio is > 1 on the first spec- 
trum of the second night and decreases until it reaches unity on 
the last spectrum (0=0.701). Significant changes are apparent 
in the Hei AAAll absorption profile, which is much broader 
and appears to be composed of two troughs, with the bluer 
trough being stronger. An additional, distinct weak blueshifted 
absorption line is also present and becomes stronger with time. 
The H6 line appears to be the result of the blending of two 
lines, showing a wide profile with two troughs. The bluer 
trough is stronger than the redder during the second night. The 
N in /14634-40-42 emission is observed as two separate peaks, 
with the redder being the strongest and the intensity of both 
increasing through the night. The He n /14686 emission feature 
appears to be symmetric compared to the first night and weak- 
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Fig. 4. B17 red spectra showing Ha variability at different 
phases taken over a period of three years. Prominent lines are 
labeled. 
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Fig. 5. Upper panel: Light curve folded on a 4.0217 day period. 
Lower panel: Equivalent widths of the He n i4686 emission 
line for the WHT data taken in 2007 August plotted on the 
same orbital period. 
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ens throughout the night. We note that (as mentioned in Section 
3) the light curve demonstrates an asymmetry - possibly due to 
mass transfer - between phases <p ~0.6-0.75. Given this is co- 
incident with the observations during night two, it is tempting 
to attribute the spectral changes (in part) to such a process. 

4.3. Spectral Classification 

Given that the current spectroscopic dataset does not span an 
entire orbital cycle we have not attempted a formal deconvo- 
lution of the blended spectra. Instead we have simply used the 
spectra closest to the primary and secondary eclipses to per- 
form a preliminary classifica tion of the components of the sys- 
tem, using the stellar atlas o f Walborn & Fitzpatrick (1990) and 
Walborn & Howarfhl d200dl) . 

The spectrum at 0=0.453 shows a He n A4541 / He i A4471 
absorption-line ratio of one, suggesting a likely spectral clas- 
sification of 07. Both Hen A4686 and Nra ^4634-40-42 are 
in strong emission which is characteristic of an Of supergiant, 
implying an initial classification of 07 Iaf. Fig. 7 upper panel 
shows this spectrum along with comparison spectra. The spec- 
trum lacksjhe Si iv /14089 line, which is unexpected for an 07 
SCifl lWarbornl tboOli mentions that a similar effect is found in 
spectra in the Small Magellanic Cloud where no Si iv /14089 is 
observed and he attributes this to the absorption and emission 
features cancelling each other. Therefore we suggest a prelim- 
inary spectral classification of 07 Iaf+ for this component. 

The spectrum at 0=0.036 (Fig. 7, lower panel) shows a He n 
/14541 / He i /14471 absorption-line ratio smaller than one, in- 
dicating a spectral type later than 07. The He n /14686 and N in 
/l/t4634-40-42 are also in strong emission, implying an Of su- 
pergiant. The spectrum also presents a stronger Si iv 4089 ab- 
sorption line compared to the spectrum at 0=0.453 although 
it is still weaker than expected for a late O star. Therefore we 
suggest a preliminary spectral classification of 09 Iaf for this 
component. 

Since the combination of spectra of the preliminary clas- 
sifications of an 07 and an 09 components would not give 
a HeI4471/HeII4541 ratio less than unity (see Section 4.2), it 
suggests there is emission in-filling of HeI4471. 

4.4. Radial Velocities 

The high resolution blue spectra show systematic night to night 
variations in the H, He and N lines, revealing significant ra- 
dial velocity shifts. Radial velocities were determined for the 
primary lines which appear not to be blended (Hen /14200, 
He ii /14541, Hy /14340) in the high resolution spectra obtained 
on the WHT observing run in 2007 August. Gaussians were 
fitted to the line profiles using the dipso emission line fitting 
command (elf; the results are discussed in Sec. 5). The lines 
for the secondary appear to be blended and were not measured. 
The radial velocities for the secondary will be measured in the 
future paper after the spectra have been dissentangled. 




Wavelength (A) 




4200 4300 4400 4500 4600 4700 



Wavelength (A) 

Fig. 7. Upper panel: B17 spectrum obtained with the WHT on 
the night of 2007 August 21, close to the secondary eclipse, 
compared with Of supe rgiant spectra from th e Digital Atlas of 
Stellar Classification d Walborn & Fitzpatricklll990l) . The spec- 
trum is most compatible with an 07Iaf classification. Double 
lines are observed even though the system is at eclipse which 
we currently cannot account for and will be further studied with 
the spectra dissentangling. Lower panel: B17 spectrum ob- 
tained with the WHT on the night of 2006 June 1 1 , close to the 
primary eclipse, compared with Of sup ergiant spectra from the 
Digit al Atlas of Stellar Classification (IWalborn & Fitzpatrick 
1990). The spectrum is more comparable to the 09 classifica- 
tion. 



5. Light and radial velocity curve analysis 

In order to reproduce the observed characteristics of the photo- 
metric light curve, w e analysed it with the 2003 version of the 
Wilson & Devinnevl (W-D; 1 1 97 lb code. For the analysis, the 
ROTSE (with an effective wavelength similar to Johnson R) 
and amateur V light curves were modeled as different datasets. 
In all cases, detailed reflection-model and proximity-effect cor- 
rections were included. Considering the spectroscopic analy- 
sis, the temperature of the primary (7? ff ) was fixed to 35 000 K, 



4 We note that while this line is not present in the spectrum of the 
07 supergiant Sanduleak 80, given the reduced metalicity appropriate 
for a SMC star it is not clear it provides a valid comparison. 



and the bolometric albedo and gravity brightening coefficients 
were set to unity, as generally found for stars with radiative 
envelopes. In addition, a circular orbit was adopted, as sug- 
gested by the equal separation between both (primary and sec- 
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ondary) eclipses, and a rotation rate synchronized with the or- 
bital period was assumed for both components. The fitting pro- 
cess was carried out iteratively until three consecutive solutions 
provided differential corrections for all the parameters smaller 
than twice their internal errors. 

Considering the characteristics in the light curves described 
above, the first runs in the modelling with W-D assumed a semi- 
detached configuration. Numerous attempts were performed 
with several mass ratios and with either the primary or the 
secondary component filling the Roche lobe. However, in all 
cases, the fits provided solutions where both stars tended to be 
in contact. Therefore, despite the properties of the lightcurve 
suggesting a semi-contact configuration we finally attempted 
an over-contact solution - the results of which are described 
here - where each component could have a different tempera- 
ture, as observed from the different depths of the eclipses. For 
this model, the mass ratio was set to q — 0.75 - as might be ex- 
pected for 07 and 09 super-giant components - noting that the 
mass ratio cannot be smaller than O.^J- The time of minimum 
(fminX the period (P), the orbital inclination (/), the temperature 
of the secondary component (7f ff ), the surface potential (£2 P ) 
and the luminosity of the primary component (L p ) were all left 
as free parameters in the fit. 

The V light curve (which has smaller photometric errors) 
clearly reveals that one quadrature is ~0.02 mag brighter than 
the other. Two main solutions were attempted to explain the ob- 
served O'Connell effect, one with an equatorial hot spot (30% 
hotter than the photosphere) on the primary component and an- 
other one with the spot being on the secondary component. The 
size and position of the spot were left as free parameters and 
also fitted at each run, although we emphasise that the size and 
temperature of the spot are strongly correlated. The fits with the 
spot on the primary component were finally adopted, since they 
provided slightly smaller errors. In addition, the position of the 
spot (oriented roughly towards the secondary component) can 
more easily be explained as the interaction of the stellar winds. 

Together with the light curve analysis, the radial velocities 
were also used to constrain the solution. Radial velocities were 
fitted separately from the light curve to avoid the larger num- 
ber of photometric observations dominate in the solution; the 
parameters fitted being the semi-major axis (a) and the sys- 
temic velocity (y). The rms of the fit is 9.5 km s _1 for the radial 
velocities, 0.023 mag for the V light curve and 0.048 mag for 
the ROTSE light curve. Unfortunately, the accuracy of the solu- 
tion obtained is strongly dependent on the mass ratio adopted. 
Nevertheless, Fig. 8 and Fig. 9 show the light curve and radial 
velocity curve with the best fit and Table[2]gives the parameters 
obtained from this solution. 

Despite our efforts we are still not completely satisfied with 
this fit. In particular, we are still unable to adequately fit the 
egress from secondary minimum nor the subsequent lightcurve 
between phases 0.6-0.9. Moreover it is expected that the tem- 
peratures of stars in over-contact binaries should be ~equal 
where a ratio of 0.85 was found between the secondary and 



Table 2. Results from the analysis of the light and radial veloc- 
ity curves: The errors shown should be considered internal er- 
rors of the fit. Any possible systematic errors are not included. 



Parameter 


Value 




T (MJD) 


4272.534 ± 0.004 




Period, P 


4.02174 ± 0.00003 days 




Inclination, i 


72 ± 1.5 




Eccentricity, e 


(Fixed) 




Mass ratio q 


0.75 (Fixed) 




Temperature ratio 


0.85 ± 0.02 




Semi-major axis (R Q ) 


50 ± 1 




Systemic velocity 


-43 ± 4 km s" 1 






Parameter 


Primary 


Secondary 


Mass (M Q ) 


60 ±5 


45 ±4 


Radius (R e ) 


22 ± 1 


19 ± 1 


log g (cgs) 


3.53 ±0.01 


3.52 ± 0.01 


RV Semi-amplitude (km s -1 ) 


257 +7 


343 ±9 


mean Teff (K) 


35000 (Fixed) 


29900 + 700 


Surface potential (Q) 


3.14 + 0.02 


3.14 (Fixed) 


M flo/ 


-9.8 ± 1 


-8.8 ± 1 



primary. As such we regard the parameters presented in Table 
2 as provisional at present. We anticipate that the determina- 
tion of radial velocities for the secondary component, which 
will directly constrain the binary mass ratio will greatly clarify 
the fundamental properties of the components and the configu- 
ration of this EB system. 

Finally, the lightcurve modelling, allows us to address the 
distance to B17, and by extention the C yg OB2 association . 
Adopting the bolometric corrections from lMartins et al.1 (|2005) 
we may use the bolometric luminosities determined above (See 
Table 2) to calculate the absolute V magnitudes of both com- 
ponents. The V band reddening w as then calculated by fol low- 
ing a similar procedure used by iNegueruela et al. 2008 (see 
Section 6) With the absolute V magnitudes, and the V band 
reddening for B17, along with the V band values for the two 
minima in the lightcurve, the distance modulus was calculated 
to be 10.9-11.3 which corresponds to a distance of about 1.5- 
1.8 kpc. While we regard these values as provisional due to the 
difficulties in modelling described above, these are consistent 
with the commonly adopted distance estim ate o f 1.7 kpc. (e.g. 



Hanson! 2003 . Torres -Dodgen et al. 1991 and Kiminki et al 



5 For a mass ratio of 0.4, the implied mass of the primary - 
considering the minimum amplitude of the radial velocity curve- is 
over 120 M Q , regardless of the light curve analysis. 



2007). Therefore they d o not agree with the distance estimate of 
900-950 pc obtained bv lLinder et ail d2009l) from the lightcurve 
modelling of Cyg OB2 #5 . 

A distance of 900 pc would imply radii that are ~half those 
currently derived from the modelling. The derived radii scale 
linearly with the semi-major axis and the semi-major axis in 
turn is strongly dependant on the assumed mass ratio. For this 
system, a semi-major axis a factor of 2 smaller would imply a 
mass ratio of 6. We consider this to be extremely unlikely, but 
not impossible. 
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Fig. 8. The V band and R band light curves of B 17 with the best 
fit model overlaid. The panels below the light curves show the 
deviations from the model. 




Fig. 9. The RV curve for the primary component obtained from 
the WHT data, with the best fit overlaid (solid line). The panel 
below shows the residuals between the observed data and the 
model. 



(2005|), and the intrinsi c (J - V)o and (K - V)q colour calibra- 
tions of lWegnerl dl994h . Using the 2MASS observed (J - K s ), 
we derived E(J - Ks). The reddening of the system was cal- 
culated using the relation Ak s - 0.61E(J - Ks), du e to the 
redde ning in the association being close to standard ([Hanson 



2003). The reddening in the V band was calcu lated using the 



relation Q.\\2A V ^A K dRieke & Lebof skv 1985). The absolute 
V band magnitude for the primary was then calculated, using 
the V band value of the secondary minimum and adopting the 
distance modulus of 1 1 .3, obt ained by averaging spectroscopic 
distances dKiminki et al.l2007l) . A semi-observational M\, a \ was 
then calculated to be -9.8 + 0.2 and -9.2 + 0.2 for the primary 
and secondary, by adding (V - K)o and the BC to the Vb> as- 
suming an uncertainty of one spectral type. With this value, the 
luminosities were calculated to be log(Li/L ) = 5.8 ±0.1 and 
log(L2/L Q ) = 5.6 + 0.1, noting that the main source of error 
is likely to be the uncertainty in the spectral type and hence 
temperature and BC. We regard these estimates as upper limits 
since it assumes that both components are completely eclipsed 
during the minima. 

The positions of the two components of B 17 in the HR di- 
agram (Fig. 10) are consistent with other known members of 
Cyg OB2 sugges ting that it too is a bona fide member, with an 
age of ~2.5Myr dNegueruela et al.ll2008l) . 

With a likely spectral type of 09 Iaf, the secondary ap- 
pears slightly more evolved than the primary and hence was 
likely the initially more massive star, with both evolving from 
very early O main sequence stars. We note that a system com- 
posed of 07 Iaf & 09 Iaf stars is also compatible with the 
other massive evolved binaries in Cyg OB2. Indeed, the 6.6 da 
Ofpe/WNL+06.5-7 binary Cyg OB2 #5 dLinder et alj|200 
appears to be remarkably similar to B17, although with the 
Ofpe/WNL star being slightly more evolved than the 09 Iaf 
secondary in B17. The slightly longer period of Cyg OB2 #5 
can be explained by one or more from (i) initial birth parame- 
ters, (ii) mass transfer from one star to the other and (iii) mass 
los t by both stars v ia stell a r wind s. 



s 



Marti ns et alJ (I2007L 120081) used observations of the 



Galactic Centre and Arches clusters to examine the evolu- 
tionary pathways of stars at galactic metallicity with masses 
in excess of 30 M o , suggesting the fo llowing evolutionary 
pathways (also see [Crowther et al ]fl995l) : 

• ~30-60M G : O Ofpe/WN9 ^± LBV -> WN8 -> WN/C 

• ~60-120M Q : O ^ Of ^ WNL + abs ^WN7 



6. Discussion 

Given the difficulties in determining a unique model fit to the 
data (Sect. 5), we regard the parameters presented in Table 2 
to be provisional. As a check on consistency we have cal- 
culated a preliminary bolometric magnitude for the compo- 
nents, following a sim ilar procedure the that employed by 
Negueruela et al.Jd 2008). Assuming that both components are 
of simi lar coloucl we adopted the 2MASS value for (J - Ks) 



used in lNegueruela et al.1 (120081) . the effective te mperature T e g 
and bolometric correction (BC) calibrations of iMartins et al 



6 The colour differences between a n Q7 and an Q 9 supergiant are 
(J - V) Q = 0.09 and (K - V) = 0.11 (Weg neill 19941) 



With an age of ~2.5 Myr and spectral classifications of 
07 and 09, B17 appears to lie at the dividing point be- 
tween the 2 evolutionary pathways. However, given that it ap- 
pears significantly less massiv e than the WN6ha + WN6ha bi 



nary W R20a (m^m^SSM jRauw etaDl2004t: iBonanos et al 



2004) we suspect that it will instead evolve to resemble Cyg 
OB2 #5 and hence to a configuration similar to GCIRS 16SW, 
composed of two coole r extreme B super giants/LBV candi- 
dates (P or b = 19.45 days; Martins et alJ2006l) . with stellar mass 
loss resulting in an eventual lengthening of the orbital period 
of B17. Indeed, assuming the system avoids merger during the 
LBV phase and, remaining bound, receives a favourable SNe 
kick to reduce the orbital separation it might briefly form a high 
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mass X-ray binary with a WR mass donor prior to the second 
SN. 

Irrespective of its ultimate fate, B17 is amongst the bright- 
est/most massive systems in Cyg OB2 and adds to the 
increasing number o f massive binaries identified within it 
( Ki minki et al. 2009, see Table 3). Similar trends for both 
a high binary fraction and the multiplicity of the bright- 
est/most ev olved cluster member s have been observed for 
Pismi s 24 dMafz ApeUanizI l2008h. NGC3603 dSchnurr et al.1 
20081) . Westerlund 1 dClark et all2008HRitchie et al.ll2009h and 
potentially the Arches ( Clark et all 2009 ). If these trends con- 
tinue it will have significant implications for the formation 
channels and relative production rates of both low and high 
mass X- ray binaries and systems c omprising of two relativistic 
objects dKobulnickv & Frverll2007h . 



-a - 



-e 



_logt-6.3 09 



logt=a.4 (rot) 



ogt=G_4 




log T c , 



Fig. 10. Updated sem i-observational HR diagram from 



Negueruela et al. (2008), based on published spectral classes 
and 2MASS JHK S photometry and a DM of 1 1 .3. The continu- 
ous l ines are non-rotating isochrones for log t = 6.2, 6.3 and 6.4 
from lSchaller et alj d 19921) and the dashed li ne is the log t = 6.4 
isochr one in the high-rotation models from lMevnet & Maedei 
(2003). The circle and the triangle show the positions of the 
primary and secondary components of B17. 



7. Summary 

Using photometric and spectroscopic data, we have demon- 
strated that B 17 is an eclipsing, double lined spectroscopic bi- 
nary comprising two supergiants with preliminary classifica- 
tions of 07Iaf and 09Iaf. The spectra are highly variable, and 
with a subset revealing features from both stars, raise the pos- 
sibility of achieving a dynamical mass determination for both 
components. Utilising both the photometric lighturve and our 



limited RV dataset we attempted to determine an initial orbital 
solution for the binary. 

Despite the morphology of the lightcurve indicating a semi- 
contact configuration we were unable to to achieve conver- 
gence for such a hypothesis and hence were forced to adopted 
an over-contact configuration. In the absence of a full RV curve 
for both system components we were forced to fix the binary 
mass ratio, and had to include the presence of a star spot to 
address the observed asymmetries in the lightcurve (which are 
likely due to the effects of binary mass transfer). However, we 
were still unable to fully fit both secondary eclipse and the 
lightcurve between orbital phase ~0.6-0.9 with such a model; 
as such we regard the modelling results presented in this work 
as provisional. We anticipate that a full RV curve for both com- 
ponents will be necessary to obtain more precise parameters 
for the system; additional data to accomplish this goal are cur- 
rently being obtained and refined analysis will be presented in 
a future paper. 

Nevertheless, the provisional distance calculation of 1.5- 
1.8 kpc obtained from the lightcurve analysis agrees with pre- 
viously published values for the distance to Cyg OB2, being 
inconsistent with t he distance of 900-950 pc determined by 
Linde r et al. (2009) with the Cyg OB2 #5 light curve analysis. 

When placed in the HR diagram, B17 appears to be con- 
sistent with the age and stellar population of the Cyg OB2 as- 
sociation. Assuming the system avoids merger, it is likely to 
evolve through an extreme B supergiant/LB V phase into a long 
period WR+WR binary configuration as mass loss via stellar 
winds increase s the orbital separation . In combination with the 
recent work of [Kiminki et al.1 (2009) and [kobulnick v & Frverj 
(2007) the results of our analysis provides additional evidence 
that Cyg OB2 has a very high fraction of massive binary stars. 
Such an observational constraint needs to be considered when 
determining the initial mass function of the association as it 
may both influence the slope of the relationship and also lead 
to a population of artificially massive stars, resulting in the in- 
flation of a putative high mass cut-off to the IMF. 
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Fig. 6. Blue spectra obtained with the WHT on the night of 2007 August 21 and 22. The main transition lines are shown. On the 
regions where the lines are blended, accuracy of the identified lines is not certain and will be clarified with spectra dissentangling. 



